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O \ Abstract 

We study the effects of noncommutativity, in the form of a Lie-algebraically deformed Poisson 
' commutation relations, on the evolution of a Bianchi type I cosmological model with a positive 

cosmological constant. The phase space variables turn out to correspond to the scale factors of 
q ■ this model in x, y and z directions. According to the conditions that the structure constants 

(deformation parameters) should satisfy, we argue that there are two types of noncommutative 
jj^ ' phase space with Lie-algebraic structure. The exact classical solutions in commutative and type I 

bJT), noncommutative cases are presented. In the framework of this type of deformed phase space, we 

investigate the possibility of building a Bianchi I model with cyclic scale factors in which the size 
^vq , of the Universe in each direction experiences an endless sequence of contractions and reexpansions. 

■ We also obtain some approximate solutions for the type II noncommutative structure by numerical 

methods and show that the cyclic behavior is repeated as well. These results are compared with 
the standard commutative case, and similarities and differences of these solutions are discussed. 

<3\ • 
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■ 1 Introduction 

'_. _ 

Since 1947 when noncommutativity between space-time coordinates was first introduced by Snyder 
[1], many efforts have been made in this area and the corresponding results have been followed by 
a number of works, the main motivations of which lie in the results of string theory [2]. Although 
the interest in noncommutative theories has been gathering pace in recent years because of strong 
motivations in the development of string and M-theories, they may also be justified in their own right 
because of the interesting predictions they have made in particle physics [3] , quantum [4] and classical 
mechanics [5] . In a general classification there are three sorts of noncommutative space-time with the 
following deformed algebra between coordinates [6] 

• Deformed Poisson algebra with canonical structure 

{%fll%v\ ^llVl (1) 

where the tensor (with constant components) 9n V is assumed to be antisymmetric. 

• Deformed Poisson algebra with Lie-algebraic structure 

{x^x u } = 6 X x x , (2) 
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where the (constant) deformation parameters 0* v are assumed to be antisymmetric to lower indices. 
• Deformed Poisson algebra with quantum space structure 

{x^,Xjy} = 9^x a xp, (3) 

where again the (constant) deformation parameters are assumed to be antisymmetric to lower 
indices. 

In general, as is clear from (l)-(3), noncommutativity in their original form imply a noncommu- 
tative underlying geometry for space-time. However, formulation of gravity in a noncommutative 
space-time is highly nonlinear, rendering the setting up of cosmological models difficult. In this 
respect, a different approach to noncommutativity is through its introduction in the phase space 
constructed by minisuperspace fields and their conjugate momenta. Alternatively, in cosmological 
systems, since the scale factors, matter fields and their conjugate momenta play the role of dynamical 
variables of the system, introducing noncommutativity in the corresponding phase space is particu- 
larly relevant [7]. This means that in the relations (l)-(3) the space-time coordinates can be replaced 
by the phase space coordinates to construct a noncommutative or deformed phase space. In this 
sense, one can introduce different kinds of noncommutativity between different dynamical variables 
of the corresponding phase space and the main aim is to study the aspects related to the application 
of noncommutativity in the framework of minisuperspace reduction of dynamics. 

In this paper we are going to investigate the impact of the deformation of Lie-algebraic type on 
the cosmological dynamics of Bianchi type I model. Since the Bianchi models have different scale 
factors in different directions, they are suitable candidates for studying noncommutative cosmology. 
Here, our aim is to introduce Lie-algebraic noncommutative scale factors in Bianchi type I space-time 
and compare and contrast their solutions to that of the conventional case. It should be emphasized 
that when we speak of noncommutativity in this work, we mean noncommutativity in the fields (scale 
factors) and not the coordinates, that is to say that we study noncommutativity within the context 
of phase space only. We shall see that the main feature of the resulting noncommutative model is 
constructing a scenario of cyclic cosmology governed by introducting a Lie-algebraic deformation in 
the phase space of Bianchi type I space-time. In such a scenario the history of the Universe (in each 
direction) is periodic, i.e., the Universe undergoes a periodic sequence of expansion and contraction 
phases. The existence of cyclic solutions in the FRW models was first introduced by Tolman in [8] 
and then the early results in this subject have been followed by a number of works, see [9] and for a 
review see [10]. In these works cyclic cosmologies which are an extension of bouncing scenarios, are 
used as an alternative for conventional inflationary cosmology. Instead of a continuously expansion 
after big-bang, the bang is replaced by a transition to an earlier phase of evolution. 

The structure of the paper is as follows: after a brief review of Bianchi type I cosmology with a 
cosmological constant in section 2, we deal with a Lie-algebraically deformed phase space in section 
3. In this section by imposing some constraints ( Jacobi identity) on the structure constants of the Lie 
algebra, we present two kinds of Lie-algebraic deformed of Poisson brackets according to which we 
write the (Poisson) commutation relations between the dynamical variables of the Bianchi I model. 
The exact solutions in type I noncommutative cases are presented. We also obtain, in this section, 
some approximate solutions for the type II noncommutative structure by numerical methods. Finally, 
in section 4 we summarize and discuss the results. 

2 Bianchi type I cosmology with cosmological constant 

In this section we make a quick review of some of the important results in the Bianchi type I model 
and obtain its Lagrangian and Hamiltonian in the ADM decomposition, for more details see [11]. 
The simplest generalization of the flat FRW cosmology is the Bianchi type I Universe which has the 
metric 

ds 2 = -N 2 {t)dt 2 + a 2 {t)dx 2 + b 2 {t)dy 2 + c 2 (t)dz 2 . (4) 
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In this metric N(i) is the lapse function and there are three functions a(t), b(t) and c(t), to be 
determined by the Einstein field equations, and are the scale factors of the corresponding Universe 
in the x, y and z directions respectively. The scale factors in different directions are allowed to vary 
independently of each other. This metric is the simplest anisotropic and homogeneous cosmological 
model which, upon making the scale factors equal, becomes the flat FRW metric. In general, the nine 
Bianchi (class A) models are the most general homogeneous cosmological solutions of the Einstein field 
equations which admit a three-dimensional isometry group, i.e. their spatially homogeneous sections 
are invariant under the action of a three-dimensional Lie group. To transform the Lagrangian of the 
dynamical system which corresponds to the Bianchi cosmologies to a more manageable form, it is 
useful to introduce the following change of variables 

a (t) = e v(t)+v(t)+V3w(t) ^ fe (f) = e u(t)+v(t)-V3w(t) ^ C (Q = e u(t)-2v(t)_ ^ 

In the Misner notation [12], the metric of the Bianchi models can be written in terms of these new 
variables as 

ds 2 = -N 2 {t)dt 2 + e 2u (V ft ^V ® co j , (6) 

where V(t) = e 3u ^ = abc is the comoving volume of the Universe and /%j determines the anisotropic 
parameters v(t) and w(t) as follows 

fiij = diag (v + y/Sw, v - V3w, -2vj . (7) 

Also, in metric (6), the one- forms uj % represent the invariant one- forms of the corresponding isometry 
group and satisfy the following Lie algebra 

dJ = -\c\ k J Nu\ (8) 

where C* fc are the structure constants. Indeed, the Bianchi models are grouped by their structure 
constants into classes A and B. Because of the difficultly in formulating the class B Bianchi models 
in the context of the ADM decomposition and canonical quantization [13], it is usually the case that 
one confines attention to the class A models where the structure constants obey the relation Cjj = 0. 
The Einstein-Hilbert action is given by (we work in units where c = h = 16nG = 1) 

S = j d 4 x^(K - A), (9) 

where g is the determinant of the metric, 1Z is the scalar curvature of the space-time metric (6) and 
A is the cosmological constant. In terms of the ADM variables, action (9) can be written as [14] 

S = J dtd 3 xC = J dtd 3 xN^h (KijK ij - K 2 + R - a) , (10) 

where Kij are the components of extrinsic curvature (second fundamental form) which represent how 
much the spatial space hij is curved in the way it sits in the space-time manifold. Also, h and R are 
the determinant and scalar curvature of the spatial geometry respectively, and K represents the 
trace of Kij . The extrinsic curvature is given by 

% = ^(% + %-^). (") 

where N^j represents the covariant derivative with respect to hij. Using (6) and (7) we obtain the 
nonvanishing components of the extrinsic curvature and its trace as follows 

' K n = -±(u + v + V3w)e 2( > u+V+V3w \ 
K22 = -jj(u + v- V3w)e 2 ^ u+v -^ w \ 

(12) 

K 33 = -^(u-2v)e 2 ^\ 
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where a dot represents differentiation with respect to t. The scaiar curvature R of a spatial hyper- 
surface is a function of v and w and can be write in terms of the structure constants as [15] 

R = C} k C l mn h a h k V + 2C) k Cih? 1 . (13) 

The Lagrangian for the Bianchi class A models may now be written by substituting the above results 
into action (10), giving 

C = ^- (-n 2 + v 2 + w 2 ) + Ne 3u (R-A). (14) 
The momenta conjugate to the dynamical variables are given by 

d£ 12 ■ 3u dC 12 ■ 3« 9C 12 . 3u 

P - = ^ = ~N Ue ^ = ^v = N Ve Pw = dv = N We (15) 

leading to the following Hamiltonian 

■H = ^Ne~ 3u (-pi +p 2 v + p 2 w ) - Ne 3u (R-A). (16) 

The preliminary set-up for writing the dynamical equations is now complete. The cosmological 
dynamics of the Bianchi models are studied in many works [11]-[16]. In the following, we shall 
consider only the simplest Bianchi class A model, namely type I. The structure constants of the 
Bianchi type I are all zero, that is C* fc = 0. It then follows from equation (13) that R = 0. Thus, 
with the choice of the cosmic time gauge N = 1, the Hamiltonian can be written as 

n = hf Zu +p " + pl ) + Ae3u - (17) 

The Poisson brackets for the phase space variables are 

{xi,Xj} = 0, {pi,Pj} = 0, {x i ,Pj} = 5 ij , (18) 

where Xi(i = 1,2,3) = v,w,u and pi(i = 1,2,3) = p v ,Pw,Pu- Therefore, the classical dynamics is 
governed by the Hamiltonian equations, that is 

f« = {^} = -T2e~ 3u p u , p u = {p u ,H} = \e~ 3u (-pI+pI+pI) - 3Ae 3u , 

v = {v,H} = ±e- 3u Pv , p v = { Pv ,H} = 0, (19) 

s w = {w,H} = j^e- 3u p w , p w = {p w ,H} = 0. 
We also have the constraint equation % = 0, from which we obtain 

e~ 3u (-Pi +P 2 v + Pi) = -24Ae 3 «. (20) 
Using this relation in eliminating p u from the first equations of the system (19) results 

il + 3u 2 = ^-, (21) 

where u> = The solution of this equation can be read as 

u(t) = \c+\]n[whhu(t-S)], (22) 

where C and 5 are integrating constants. Since 5 only corresponds to a shift in the big-bang singularity, 
we can put it equal to zero. The last two equations of the system (19) indicate that v(t) and w(t) 
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obey the equations of motion v = ^Po v e~ 3u and w = j2Po w e~ 3u , where we take p v = po v = cons, 
and p w = po w = cons. Upon substituting the relation (22) into these equations, they can easily be 
integrated to yield 



17(f) 



POv 

12e c co 



In 



tanh — 
2 



w(t) 



POv 



l2e c uo 



In 



tanh — 
2 



(23) 



Now, these solutions must satisfy the zero energy condition, % = 0. Thus, substitution of equations 
(22) and (23) into (20) gives a relation between integration constants as 



C = ln 



POv + POw 



(24) 



From these relations one can calculate all of the physical quantities of observational interest in cos- 
mology. In analogy with a FRW Universe, we can define a Hubble parameter corresponding to each 
of the three different directions 



H r 



7 1 ~ ± C ) 

a b c 



(25) 



in terms of which the average Hubble parameter is defined as 



H = -(H a + H b + H c 



u. 



(26) 



The directional Hubble parameter Hi measures the expansion rate of the Universe in the direction Xi 
while the average Hubble parameter H measures its volumetric expansion rate. We will also define the 
expansion scalar B = 3H, the shear scalar a 2 = ^ (j2i=i Hf — 3H 2 ^j = 3 (v 2 + w 2 ), which measures 

the degree of anisotropy of the space-time and the average deceleration parameter q = ^(-jj) — 1, 
where as is well known is indicated by how much the expansion of the Universe is slowing down. If 
the expansion is speeding up, for which there appears to be some recent evidence, then this parameter 
will be negative. Now, let us return to the variables a,i{t) \a\{t) = a(t), ci2(t) = b(t) and a^t) = c(t)] 
using the transformation (5), in terms of which we obtain the corresponding scale factors for the 
Bianchi I cosmology as 



ai(t) 



\/P0v +Plv 

2uj 



1/3 



smh» + 1 — cosh a Sl — , 
2 2 ' 



(27) 



where the constants s,; are defined as 



POv + V^POw POv ~ V^POw 
Sl = = , S 2 = ; S 3 



3\/Pov +P0u 

and satisfy the following relations 



SJPov +pL 



-2pov 



3JP0v + POn 



(28) 



si + s 2 + s 3 = 0, sl + sl + s§ 



(29) 



From these equations we obtain, the directional Hubble parameters 



<H SiU w 

H i = — = — 1 cothwf, 

ctj smh cot 3 



the comoving volume 



(30) 



(31) 
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Figure 1: Top: Behavior of the scale factors a(t), b(t) and c(t) (solid lines) and their corresponding Hubble parameters 
(dashed lines) as a function of time. Bottom: The comoving volume V(t) = e 3u — abc (solid line) and its corresponding 
Hubble parameter (dashed line), the shear scalar and the deceleration parameter. The figures are plotted for numerical 
values A = 1, po v = 1 and po w = 2. 



the average Hubble parameter 
the shear scalar 

and the deceleration parameter 



H(t) = ^cothwi, (32) 

■*> - xisL? < 33 > 

«(') = — ~ L < 34 > 

cosh cut 

In figure 1, we have shown the behavior of the scale factors, Hubble parameters, comoving volume, 
shear scalar and deceleration parameter for typical values of the parameters. It is seen that the 
evolution of the Universe begins with a big-bang singularity at t = with a high degree of anisotropy 
and then follows an expansion phase in all directions. Since we have negative acceleration (positive 
deceleration parameter) at early times, the Universe decelerates its volumetric expansion in this era. 
On the other hand, the positive acceleration (negative deceleration parameter) is occurred for late 
times which means that the Universe currently accelerates its volumetric expansion. While the volume 
of the Universe and its scale factors in different directions increase monotonically, the shear scalar 
decreases monotonically and tends to zero as t — > oo. This means at late time the model approaches 
that of the flat FRW Universe with cosmological constant, i.e., the Universe eventually evolves to a 
phase where it is close to a de Sitter solution. 



3 Lie-algebraically deformed phase space 

We now focus attention on the study of the effects of Lie-algebraic noncommutativity on the cos- 
mological dynamic of the Bianchi type I model described above. This kind of noncommutativity is 
described by a deformed Poisson algebra between the phase space variables with the Lie-algebraic 
structure as follows 

{xi,Xj} = OijXk, (35) 

where the constants 9fj are assumed to be antisymmetric with respect to their lower indices. To 
guarantee the Jacobi identity 



{{Xj, Xj}, Xfc} + {{Xj, Xfc}, Xj} + {{Xfc, Xj}, Xj} — 0, 



(36) 



these parameters are also subject to the following conditions 

al am , al am , al am n (oi7\ 

VijVlk + Vkfllj + VjkVii - U. (67) 
Following [17], we consider the Poisson brackets between the dynamical variables and their momenta 



as 



{xi,pj} = Sij + 0\jX k + OijPk, {pi,Pj} = 0, (38) 



in which the new noncommutative parameters 9 and 9 with vanishing diagonal elements are introduced 
and furthermore, we assume that the momenta (Poisson) commute with each other. For our model 
at hand, Xi(i = 1,2,3) = v,w,u and pi(i = 1,2,3) = p v ,p w ,p u , as before. In the phase space with 
coordinates ff = (x, p), in addition of the identity (36), the following Jacobi identities should also 
hold 

{{xi,Xj},p k } + {{xj,p k },x,i} + {{pk,Xi},Xj} = 0, (39) 



{{xi,Pj},p k } + {{Pj,Pk},Xi} + {{p k ,Xi},pj} = 0. 
These relations impose some more constraints on the noncommutative parameters as 

0, 



(40) 



*S + Oik ~ G) k 



qI am I al am 
hflk + V jk V H 



al am al am 



amal al am , al am n 



nl am 
y jk°il 



(41) 



ql am 

hfik 



ql am 
j ik°lj 



o, 

0. 



Now, the equations (37) and (41) should be solved to determine the noncommutative parameters. 
Here we will not deal with the details of the methods through which one can find these solutions and 
only refer to the Refs. [17] where the following two kinds of solutions to the equations (37) and (41) 
are proposed. 

• Type I: For this type of solution, the nonvanishing noncommutative parameters are as follows 
(the antisymmetric counterparts of Ofj should be also considered) 



e 2 



13 



9 1 

'23 



9 1 

7 32 



52 
'31 



(42) 



In this case the Poisson brackets of the phase space variables correspond to our cosmological setting 
read as 

{v,w} = 0, {w,u} = —6v, {u,v} = —Ow, 
{u,p u } = 1, {u,p v } = -9p w , {u,p w } = 6p v , 

(43) 

{v,p u } = 0, {v,p v } = l, {v,p w } = 0, 

{w,p u } = 0, {w,p v } = 0, {w,p w } = 1 

• Type II: For this type of solution, we have the following nonvanishing noncommutative param- 
eters (and also the antisymmetric counterparts of 9, 



9 2 



13 



9 1 



9 1 

7 32 



ij J 



52 
'31 



9 2 



9 1 



(44) 
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In this case the Poisson brackets of the phase space variables correspond to our cosmological setting 
read as 

{v,w} = 0, {w,u} = —9v, {u,v} = —9w, 
{u,p u } = 1, {u,p v } = -9w - 9p w , {u,p w } = -9v + 6p v , 
{v,p u } = 0, {v,p v } = l, {v,p w } = 0, 



(45) 



{w,p u } = 0, {w,p v } = 0, {w,p w } = 1 

The Poisson bracket for any phase space function can be straightforward and obtained from (43) and 
(45) and reads 

„ f) f f)n 

(46) 



{ / (x,p), 5 (x,p)}={^,^}|^^ ! 



drj A dr\ l 

where fj = (x, p) = (v, w, u,p v ,p w ,p u ). The time evolution of such a function is thus given by 

df 



dt 



{/,«}■ 



(47) 



For deformed phase space of type I the cosmological dynamics of the Bianchi-I model are 
ii = {u, %} 



1 & 3 u ~. 

12*- JrU 5 



v = { V ,H} = 9w [-|e- 3 « (-pl+pl + p 2 w ) + 3Ae 3 " 



+ 12 e ^ U Pvi 



yj = {w,H} = -ev [- |e- 3u (-pi +p 2 v + pi) + 3Ae 3 " 



{p u ,H} = -3Ae 3w + le~ 3u {-p 2 u +p 2 v + p 2 w ) , 



Pv = {Pv,H} = 9p w f-k" 3 " {-pi + p 2 v + pi) + 3Ae 3 « 



12^ Pwi 



(48) 



p w = {p w ,n} = -9p v [-|e" 3u {-Pl+Pl +Pl) + 3Ae 3u 

As before, assuming the full Einstein field equations hold, this implies that the corresponding Hamil- 
tonian must vanish, that is % = 0, which yields 



\e- Zu (-pl+pl+pl) =-3Ae 3 «. 



(49) 



We see that the noncommutativity does not affect the dynamics of u and p u . Therefore, we have 



u(t) = \-C + - ln[sinhwt], p u (t) = —4:0je c coshwt, 



(50) 



which means that the volumetric expansion of the Universe given by V{t) = e 3 "W is the same as the 
ordinary Bianchi I model described in the previous section. Substituting these results into the rest 
of the equations of the system (48), we are led to 



v = 6A9e c w sinhwt + 



12e a sinhu^ 1 " 



w = — 6A9e c v sinh cot + , c 1 



12e v sinhwt^" 



p v = 6A9e c p w sinhcji, 
Pw = — 6A9e c p v s'mhuit. 



(51) 
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From the third and fourth equations of this system, we obtain 

p v - ujp v coth ojt + 36A 2 6 2 e 2C p v sinh 2 ut = 0, (52) 
where after integration we get the following expressions for p v and p w 

( 6A ^ eC , \ ( GAee ° i \ 
Pv(t) = POv cos I cosh out I + pq w sin I cosh ut I , (53 J 

/6A0e c , \ f6A0e c , \ 
Pw(^) = — Pov sm coshwt +po«icos coshwt , (54) 

V w / V w / 

in which we have set the integration constants as pq v and pq w to achieve the correct results in the 
limit 6 — > 0. With these results at hand we can integrate the two first equations of the system (51) 
to obtain 

( 6Aee ° , \ ^ ( QAOe ° , \ 
v(i) = Ci cos [ coshwt I + 62 sm I coshwt I + 



1 / , uft\ [ (QMe c \ . (QK6e c , \1 

m I tanh — J < po^ cos I cosh I + pow sm I cosh ujt I > , (55) 



12eC 



u;(t) = C2 cos coshwt — Ci sm coshwt + 

V w / V w / 

' / , wt\ f /6A^e c , \ . (QK9e c , \1 

m ( tanh — ) pow cos I cosh ut I — po*j sm I cosh cot I > , (56) 



12e c w V 2 

where Ci and C2 are integration constants. The requirement that the deformed Hamiltonian con- 
straints (49) should hold during the evolution of the system leads again to the relation (24) between 
integrating constants. In figure 2 we have plotted the cosmological functions corresponding to the 
above solutions. As is clear from this figure the volumetric expansion, average Hubble and deceler- 
ation parameters of the deformed model are exactly the same as the usual Bianchi I model. This is 
because these parameters come from the function u(t) which is not affected by the deformation on the 
phase space. However, while like the usual model, the volume of the Universe increases continuously, 
the evolution of the scale factors in three directions is quite different in comparison with the usual 
Bianchi I model. We see that the scale factors behave cyclically, i.e., they undergo a periodic sequence 
of contractions and expansions phases. In contrast to the conventional big-bang inflationary cosmol- 
ogy described in the previous section in which the Universe begins with the big-bang singularity and 
expands forever, here we have an alternative in which the bang is replaced by a transition to an 
earlier phase of evolution. In each cycle of such a cosmology the scale factors begin their evolution 
with a big-bang and end in a big-crunch, only to emerge in a big-bang once again. Cyclic cosmologies 
have been extensively studied in literature in an attempt to solve some problems of the standard 
cosmological model such as the flatness problem, the horizon problem, the coincidence problem, the 
problem of initial conditions etc. [10]. In general, a cyclic Universe can be generated by adding a 
matter or a field that will produce a bounce, and then explore what conditions are to be imposed on 
it to produce oscillations. In the above analysis we provided an alternative method to construct a 
Bianchi type I cosmology with cyclic scale factors by introducing a deformation of Lie-algebraic type 
on the phase space of this model. In our model although the expansion of the whole volume of the 
Universe is speeding up, for which there appears to be some recent evidence, from the evolution of 
the scale factors point of view, the history of the Universe is periodic and all of the key events that 
construct the large scale structure of the observable Universe occurred a cycle ago. An important 
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Figure 2: Top: Behavior of the scale factors a(t), b(t) and c(t) (solid lines) and their corresponding Hubble parameters 
(dashed lines) as a function of time. Bottom: The comoving volume V(t) = e u — abc (solid line) and its corresponding 
Hubble parameter (dashed line), the shear scalar and the deceleration parameter. The figures are plotted for numerical 
values A = 1, po v = 1, Pow = 2, = 0.25 and Ci = C*2 = 1. 



feature of this model is that its shear scalar which measures the degree of anisotropy in the space-time, 
behaves nonmonotonically. For small values of t, it decreases and tends to zero and then it increases 
monotonically, eventually diverging as t — > oo. Therefore, the model under consideration does not 
satisfy the conditions for approaching isotropy. Since this behavior is not consistent with late time 
observations, it seems an additional mechanism is needed to modify this feature of the model. This 
is the subject of our forthcoming work and thus we do not deal with this issue here. 

For deformed phase space of type II, using the (Poisson) commutation relations (45), the cosmo- 
logical dynamics of the Bianchi-I model after applying the Hamiltonian constraint are 

u = {u,n} = - j^9e~ 3u (wp v + vp w ) - ^e" 3 >„, 
v = {v,H} = QA9we 3u + ^e- 3u Pv , 
w = {w,H} = -GA6ve 3u + ^e" 3 > w , 

(57) 

Pu = {Pu,n} = -6Ae 3 «, 
Pv = {Pv,n} = 6A(9w + 6 Pw )e 3u , 
^ p w = 6A(9v -9p v )e 3u . 

We see that this kind of deformed cosmology forms a system of nonlinear coupled differential equations 
which unfortunately cannot be solved analytically. In figure 3, employing numerical methods, we have 
shown the approximate behavior of the corresponding cosmological functions for typical values of the 
parameters and initial conditions respectively. From this figure it is seen that the cyclic behavior 
of the corresponding cosmology is almost like the case of type I deformation except that here the 
average deceleration parameter also behaves oscillatory which means that the Universe experiences a 
periodic sequence in which it decelerates and accelerates its expansion alternatively. 

4 Conclusions 

In this paper we have studied the possibility of emergence of a Bianchi type I cyclic cosmology in the 
framework of a deformed phase space of Lie-algebraic type. In the proposed model the phase space 
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Figure 3: Up: Qualitative behavior of the scale factors a(t), b(t) and c(t) (solid lines) and their corresponding Hubble 
parameters (dashed lines) as a function of time. Bottom: The comoving volume V(i) = e 3u — abc , the shear scalar 
and the deceleration parameter. The figures are plotted for numerical values A = 1, 6 = 0.5 and 9 = 0.5. We take the 
initial conditions w(0) = v(0) = w(0) = —1, po u = —2, po v — 1.7 and po w = 1. After examining other sets of initial 
conditions, we verify that this behavior repeats itself. 

coordinates obey the (Poisson) commutation relations with a Lie-algebraic structure, the structure 
constants of which play the role of the deformation parameters. These constants are subject to con- 
ditions to guarantee the Jacobi identity which by solving them the noncommutative (deformation) 
parameters can be found. In this line we have followed Rcfs. [17] that offer two types of noncom- 
mutative parameters which satisfy the corresponding constraints. We investigated the Bianchi type I 
cosmology based on these two kinds of deformed phase space and compare the results with the con- 
ventional model. For the deformed model of the first kind we showed that while the behavior of the 
whole volume of the Universe and its deceleration parameter are the same as the consensus Bianchi 
I model, the evolution of its scale factors is quite different. In spite of the usual Bianchi I model in 
which the three scale factors increase monotonically, in the Lie-algebraically deformed model they 
evolve cyclically. This means that in the scenario proposed by the deformed model the scale factors 
of the Universe undergo an endless sequence of epochs which begin with a big-bang and end in a 
big-crunch, i.e., the size of the Universe in each direction bounces from contraction to reexpansion 
alternatively. We saw that this behavior causes a high degree of anisotropy as t — > oo in contrast to 
the late time observations and therefore it seems an additional mechanism is needed to modify this 
feature of the model. For the second kind of deformation parameters, since we could not solve the 
system of field equations analytically, we investigated the subject by employing numerical methods. 
Our numerical analysis in this case showed again an oscillatory behavior not only for the scale factors 
but also for the total deceleration parameter. 
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